Comparative study of tight-binding and ab initio electronic structure calculations 
focused on magnetic anisotropy in ordered CoPt alloy 
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An empirical multiorbital (spd) tight binding (TB) model including magnetism and spin-orbit 
coupling is applied to calculations of magnetic anisotropy energy (MAE) in CoPt Llo structure. A 
realistic Slater-Koster parametrisation for single-element transition metals is adapted for the ordered 
binary alloy. Spin magnetic moment and density of states are calculated using a full-potential 
linearized augmented plane-wave (LAPW) ab initio method and our TB code with different variants 
of the interatomic parameters. Detailed mutual comparison of this data allows for determination 
of a subset of the compound TB parameters tuning of which improves the agreement of the TB 
and LAPW results. MAE calculated as a function of band filling using the refined parameters is 
in broad agreement with ab initio data for all valence states and in quantitative agreement with 
ab initio and experimental data for the natural band filling. Our work provides a practical basis 
for further studies of relativistic magnetotransport anisotropies by means of local Green's function 
formalism which is directly compatible with our TB approach. 



I. INTRODUCTION 

Ordered CoPt alloys have been studied widely as they 
hold potential for applications in high density mag- 
netic recording due to the combination of exchange 
and spin-orbit interactions giving rise to large magnetic 
anisotropies. Tunnelling magnetoresistance,^ tunnelling 
anisotropic magnetoresistance^ or spin pumping^ have 
been demonstrated in CoPt based devices. Our general 
objective is to develop an efficient numerical model al- 
lowing us to study ground state and in future studies 
also transport properties of spintronic devices based on 
CoPt or other intermetallic compound with large mag- 
netic anisotropy energy (MAE). 

MAE has been calculated in CoPt Llo structure us- 
ing ab inito methods!^ The stability of various bulk 
CoPt structures^ has been studied recently also by ab 
inito methods. On the other hand, quantum transport 
in micro-devices is typically described within the Green's 
function formalism assuming an expansion of the elec- 
tronic states on a local basis set which facilitates parti- 
tioning of the system. The tight binding (TB) descrip- 
tion of electronic structure provides a good foundation 
for subsequent simulations of magnetoresistance in tun- 
nelling or ohmic regime as it assumes a local basis set!^ 

TB schemes applied in modelling of magnetotrasport 
phenomena range from empirical or semi-empirical 
(charge self-consistent) models^Q] to tight binding lin- 
earized muffin-tin orbital (LMTO) modet^H^ combin- 
ing density functional theory (DFT) with the TB ap- 
proach. The state-of-the-art magnetic DFT-based TB 
schemeii^Eni a ll ow for simulations of complex systems 
(impurities, structural defects, surfaces) that are beyond 
practical capability of ab initio calculations. In this 
work, we employ an empirical two-center Slater-Koster 
TB modeP^ following the Harrison approach^ recently 



further developed by Shi and Papaconstantopoulos^ to 
investigate MAE of a bulk ordered transition metal al- 
loy and compare the results to density functional theory 
(DFT) calculations performed using the full-potential lin- 
earized augmented plane-wave (LAPW) program pack- 
age WIEN2k. Thereby we test the transferability of the 
TB parameter:pni obtained by fitting to ab initio band 
structures of single-element solids of the two atoms form- 
ing our ordered alloy. Rather than finding a new full 
set of TB parameters by fitting to the DFT band struc- 
ture of the compound, we determine a minimal subset 
of interatomic parameters which influence the spin mag- 
netic moments and projected density of states (DOS) in 
a transparent way and tune these parameters to improve 
the agreement of TB and DFT results. MAE calculated 
in a narrow range of refined values of the interatomic 
parameters is in quantitative agreement with DFT. We 
believe that our model and parameters provide good ba- 
sis for future simulations of magnetotransport in CoPt 
and other multilayer structures. 

Our paper is organized as follows: The TB model in- 
cluding the exchange and spin-orbit interaction is intro- 
duced in Sec.|nJ TB parameters according to Shi and Pa- 
pacons tant opoulos with our modifications are discussed 
in Sec. Ill First comparison of experimental MAEpS ab 
initio, and TB results is presented in Sec. IV Sec. [V] de- 



scribes the refinement of TB parameters towards agree- 
ment of spin magnetic moment and DOS with ab initio 
results, it also presents corresponding MAE in compari- 



son with DFT; Our work is briefly summarised in Sec. VI 



II. TIGHT BINDING MODEL 

Our TB model is a variant of the Linear Combination 
of Atomic Orbitals method for a periodic crystal where 
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the basis has the form of Bloch sums of atomic like wave- 
functions: 



(rlaak) = 



1 



N-l 



where k is the Bloch wave vector, a is the atom index 
in the unit cell, a is the atomic orbital quantum num- 
ber, N is the number of unit cells (or atoms if there is 
only one atom per unit cell), n is the unit cell index, 
R n is the unit cell vector, and p a is the position vec- 
tor of the atom a in the unit cell. We consider an or- 
thonormal set of wavefunctions 4> aa (r) constructed from 
the atomic orbitals (with the angular part expressed 
in terms of cubic harmonics given in Eq. (All) using 
Lowdin's orthonormalization procedure so the overlap 
matrix (b(3k\aak) = S^p^a- Our single-electron Hamil- 
tonian has three components as follows: 



H — H band + H stoner + Hgo ■ 



(2) 



The first component Hband contains the kinetic energy 
and a superposition of atomic potentials centered at each 
site: 



V(r) 



V a at (r-R n - Pa ) 



(3) 



where the indices n and a run through all unit cells and 
all atoms in a unit cell, respectively. The potentials 
V a at (r) centred at each site are spherically symmetric so 
the wavefunctions 4> aa can be specified by the usual an- 
gular momentum quantum numbers. 

The matrix elements of the non-magnetic Hamiltonian 
term H ba nd can be written in terms of on-site energies e aa 
and hopping integrals Ebp >a a(Pn) which depend, in the 
two-center approximation proposed by Slater and Koster, 
only on the intersite position vector p n = R n + p a — Pb 
as follows: 



'bp, act 



(bf3k\H b and\aak) = 5 b p iaa e ao 

■{Pn) = J ' dr<j)* bfj {r)H n (r)<p aa (r - p n ), (4) 

where the Hamiltonian contribution H n (r) accounts for 
interactions between atomic sites a and b in the same 
(Jin = 0) or neighbouring unit cells (R n ^ 0). 

Following the TB scheme, we parametrize the Hamil- 
tonian matrix instead of performing the integration of 
Eq. Q. The hopping integrals Ebp^ aa (p n ) can be ex- 
pressed in terms of Slater-Koster parameters (V ssa , V spcr , 
V s da, V ppa , Vpp-x, Vpda, V p dn, Vddv , Vddir , VddS , etc.) ob- 
tained usually by fitting to ab initio calculations. See 
Eqs. (Bl) and (B2) in the Appendix for the so called 



Slater-Koster tables listing the s, p, and d hopping inte- 
grals used in this work. 

The on-site energies e aa could be approximated by 
the atomic values as in the original Harrison approach. 
Instead, we follow Shi and Papaconstantopoulo:p3 who 



keep the on-site fc-independent matrix diagonal but use 
e aa obtained by fitting to ab initio data. Hence, both the 
on-site energies and Slater-Koster parameters are inputs 
of the model as we discuss in more detail in Sec. IIHI 

Since the Hamiltonian has the same periodicity as the 
basis functions in Eq. ([!]) it is diagonal in the fc-vector. 
The dimension of our Hamiltonian matrix is given by 
the number of atoms in the unit cell and valence orbitals 
considered for each atom. We use s, p, and d orbitals 
to model CoPt and there are only two atoms in the unit 
cell of the Ll structure. The sum in Eq. Q runs over 
a limited set of neighbouring sites due to localization of 
functions <p aa around site a. 

The Hamiltonian described so far is non-magnetic. 
Now we double our Hilbert space by including the 
spin degree of freedom \aak) — > \aa^k) and add a k- 
dependent term Hstoner to account for the ferromag- 
nctism in our system: 



o it < r 

\aa£k) — -Sbp^aalaafn ■ 



(5) 



where I aa and Ib/3,aa & re the on-site and hopping Stoncr 
parameters, respectively, rh is the magnetization unit 
vector constant at all sites, and cr is the vector of Pauli 
matrices. We derive our Stoner parameters from the 
exchange-split on-site energies and hopping integrals: 



-a«4-5 -*-b{3,aaL 



E, 



b/3t,aat 



E, 



ing on-site and Slater-Koster parameters fitted indepen- 
dently for the spin-up and spin-down states in the ab- 
sence of spin-orbit coupling PS Using the averages of spin- 
up/down parameters (non- magnetic structure with e aa 
and Ebp jaa ) and their differences (Stoner parameters) al- 
lows us to rotate the magnetization direction according 
to Eq. ([5]) when the spin-orbit coupling is considered. 

In order to account for the MAE we add the spin-orbit 
coupling term Hso m its atomic fe-independent form to 
the on-site terms of our Hamiltonian: 

(b(3(k\H so \aa£,k) = 5 b , a Kp (6) 

where X a p,aa is the spin-orbit parameter for orbitals a, 
/3 (p does not couple with d) and site a, L is the orbital 
momentum operator, and S is the spin operator. The 
matrix elements of Hso m the basis of cubic harmonics 



are given in Eqs. (A2) and (A3) in the Appendix 



Finally, we do not introduce explicitly any Hamilto- 
nian term controlling the charge transfer between sites 
occupied by different atoms. Prior to our calculation we 
shift all on-site energies of Co with respect to Pt so that 
the Fermi energies calculated for pure Co and pure Pt 
are equal. We check that the local charges on Pt sites 
and on Co sites in CoPt Llo structure are in agreement 
with the LAPW results within the error-bar caused by 
charge located outside of atomic spheres used by the ab 
initio method. 
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III. PARAMETRIZATION 

As mentioned above, our model relies on input param- 
eters that are obtained by fitting band structures and 
total energies to ab initio results. Extensive parameter 
sets^ are available for bulk single-element metals. They 
assume non-orthogonal basis set, interaction to higher or- 
der neighbours, and reproduce ab initio data with great 
accuracy. Our aim is to study TB models of more com- 
plex systems such as ordered binary alloys (this work) 
and heterostructures (future work) suitable for explor- 
ing the physics of relativistic equilibrium and potentially 
also transport phenomena in these systems. Therefore we 
prefer smaller, more transferable sets of parameters as- 
suming interactions only up to third nearest neighbours. 

We use a parametrisation by Shi and Papaconstan- 
topouloiP^ which further develops the Harrison ap- 
proach.^ Harrison expressed the two-centre Slater- 
Koster parameters V a p 1 {p) as functions of the inter- 
atomic distance p — \p\, an effective radius of the d or- 
bital Td which is characteristic to each transition metal, 
and constants riapy which are universal for all elements 
and lattice structures: 

V a p 7 (p) = 

mp z 

h 2 r 3 J 2 

Vady (P) — Vad-y 7 / 9 ; 

mp'' z 
h 2 r 3 

Vddry(p) = Vddy f, (7) 

mp° 

where a and /? run through the orbitals s and p. Values of 
VaPy are listed in Ref. |23j . In case of transition metals, 
Harrison used only the s and d orbitals so there were 
only two parameters specific to each element, the d-band 
width given by and the on-site energy of the d orbitals 
with respect to the s orbitals. 

Shi and Papaconstantopoulos significantly improved 
the ability of the Harrison parametrisation to produce 
accurate numerical results for the band structure while 
keeping the form and universality of the Slater-Koster 
parameters given in Eq. Q. This is accomplished by: 
1) Replacing the atomic energies by on-site energies fit- 
ted to Augmented Plain Wave (APW) calculations; 2) In- 
cluding the p orbitals into the basis set; 3) Modifying of 
the sp Slater-Koster parameters by introducing a dimen- 
sionlcss parameter 7 S as follows: 

Va/3 7 (p) = VaM 2 ' ( 8 ) 

mp z 

4) Obtaining new prefactors rj a pj and radii by simul- 
taneously fitting the APW energy bands of 12 transition 
metals at the equilibrium lattice constants of the partic- 
ular element. The new parameters reproduced APW en- 
ergy bands and density of states (DOS) remarkably well, 
not only for the 12 elements originally fitted, but also for 
the rest of the transition metals, the alkaline earth and 



the noble metals as shown in Ref. [53]. This parametri- 
sation assumes an orthogonal basis set and interaction 
to second (fee) or third (bee) nearest neighbours. The 
on-site and Slater-Koster parameters are exchange-split 
in case of the ferromagnetic metals. 

In this work we build on the results of Shi and Papacon- 
stantopoulos and test the transferability of their TB pa- 
rameters to an ordered binary alloy. We use the param- 
eters for single-element bulk metals^ with the following 
modifications: 1) The interatomic Slater-Koster parame- 
ters between Co and Pt atoms are set to a geometric av- 
erage of the elemental values: V^ Pt = JV^^ fol- 
lowing the work of Ballhausen and GraypS Using the ge- 
ometric average is also in line with the LMTCP3 method; 
2) The exchange-splittings of the on-site and Slater- 
Koster parameters enter our model through the on-site 
(Iaa) and hopping (I a a,bp) Stoner parameters, respec- 
tively. We increase Ico,d and introduce non-zero Ipt,d to 
reproduce the LAPW spin magnetic moment and DOS of 
CoPt more accurately; 3) We vary the on-site energy of 
unoccupied p orbitals (ept, p ), which extend the original 
Harrison's set of parameters and are specific to a partic- 
ular single-element crystal, in order to further improve 
the agreement with spin magnetic moment and DOS ob- 
tained by LAPW and explore the dependence of MAE in 
the CoPt compound on this parameter. 

Finally, we add atomic spin-orbit coupling parameters 
obtained by numerical Hartree-Fock calculations based 
on the Dirac-Coulomb HamiltoniarP^to the Pt sites. The 
5g? orbital has an atomic value Xpt,d = 0.0445 Ry and the 
3c? orbital of Co has Xco,d = 0.0063 Ry. We neglect the 
spin-orbit coupling in 6p-orbitals of Pt and 4p-orbitals of 
Co to keep the number of input parameters low. 



IV. COMPARISON OF TIGHT BINDING AND 
AB INITIO RESULTS 

We check the validity of our extensions of the Harrison 
TB model and the transferability of Shi and Papacon- 
stantopoulos' parametrization to CoPt ordered alloy by 
comparing the band structure, spin magnetic moment, 
DOS, and MAE calculated using our TB code and a well 
established ab initio co de, th e full-potential relativistic 
LAPW package Wien2KPESI 

In CoPt LTo structure, fee lattice sites are occupied 
by alternating layers of Co and Pt atoms and the lattice 
constant perpendicular to layers c is smaller than the 
in-plane lattice constant a. Throughout this work we 
use the experimental lattice constants: a = 7.19 a.u, 
c = 7.01 a.uP 

The Slater-Koster parameters for single-element fee 
crystals^ were fitted assuming interaction to second 
nearest neighbours. However, we intend to account for 
MAE caused by hybridisation of magnetic 3c? orbitals on 
Co and spin-orbit coupled 5d orbitals on Pt. The cou- 
pling between the Co and Pt sites in the LTo lattice in 
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(a) Second nearest neighbours 
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(b) Third nearest neighbours 

FIG. 1. (Color online) Comparison of TB (continuous red 
lines) and LAPW-LSDA (dashed black lines) band structures 
of CoPt LiO structure with magnetization along the [001] 
axis. Second (a) and third (b) nearest neighbours together 
with "atomic parameters" are assumed in TB whereas LAPW 
bands are the same in both plots. 



our two-center approximation is formed by 8 (16) second 
(third) nearest neighbour hopping integrals. Therefore, 
to enhance the ability of our model to capture the Co-Pt 
hybridisation, we add the third nearest neighbours to the 
sum in Eq. Q. 

Fig. [T] shows band structures calculated using our TB 
model assuming second (a) and third (b) nearest neigh- 
bours compared to LAPW results with local spin den- 
sity approximation (LSDA) of the exchange-correlation 
potential. At this stage we do not modify the TB pa- 
rameters derived for single- element metals except taking 
geometric average of Slater-Koster parameters between 
Co and Pt. We denote this default set of TB parameters 
as the "atomic parameters" . 

There is a good agreement of the bands especially in 
the vicinity of the Fermi energy. TB bands assuming 
second nearest neighbours match LAPW bands also in 



the region of s and p states. However, the bottom of the 
s-band is very far from the Fermi energy and the p-band 
has only an auxiliary role in our model so we can conclude 
that summation to both second and third neighbours has 
potential for successful simulation of the MAE in CoPt 
LIq structure. 

Now we turn our attention to the spin magnetic mo- 
ment, spin-resolved DOS and MAE calculated using sec- 
ond nearest neighbours together with the "atomic pa- 
rameters" and compare these TB results to LAPW with 
LSDA and spin-polarised generalized gradient approx- 
imation (GGA). It is known that GGA and LDA are 
suitable for 3d and 5d transition metals, respectively, so 
it is not clear which approximation is more appropriate 
for the CoPt compound. Therefore, we try to view the 
difference between TB and LSDA results in the context 
of the difference between calculations using LSDA and 
GGA approximations. 

Fig. [2]ja) shows the spin magnetic moment on Co and 
Pt sites per formula unit (f.u.) for energies ranging 
through the whole valence band for magnetization along 
the [001] axis. (Our [001] axis is set perpendicular to 
the alternating Co and Pt atomic planes and the nearest 
in-plane neighbour lies on the [110] axis.) The Wien2K 
code places all atoms of the unit cell in non-overlapping 
spheres leaving some charge in the interstitial region. We 
set the radii of the atomic spheres to rco = 2.2 a.u. and 
rp t = 2.4 a.u. for all LAPW calculations causing our 
spin magnetic moment in the interstitial region to be less 
than 5% of the total spin magnetic moment. GGA and 
LSDA give very similar values in case of Pt and deviate 
slightly in case of Co in the middle of the (i-band. TB 
spin magnetic moment on Pt is lower than the LAPW 
prediction in the whole energy range, whereas TB values 
for Co match the GGA very well except around the Fermi 
level where Co is predicted to be less spin-polarized by 
TB than by LSDA and GGA. The total spin magnetic 
moment is 2.27 /is/f.u. using GGA and 2.21 /ie/f.u. 
using LSDA which is in good agreement with the experi- 
mental value 2.4 /zs/f.u. TB prediction of the total spin 
magnetic moment is 1.92 /x#/f.u. Very similar results are 
obtained for magnetization along the [110] axis so we do 
not plot them. 

Fig.^b) shows DOS projected on spin-up/down states 
again for magnetization along the [001] axis and summa- 
tion to second nearest neighbours. LSDA and GGA are 
in excellent agreement in case of spin-down but differ 
slightly in case of spin-up. The DOS calculated by TB 
has unexpected peaks close to the bottom of the d-band. 
A separate calculation of the Pt-component of DOS and 
the fact that these peaks for spin- up/down are not mu- 
tually shifted in energy suggests that they correspond to 
Pt-states. Wc hypothesize at this stage that an increased 
coupling between Co and Pt sites could remove these un- 
realistic peaks. In general, the spin-up part of the valence 
band calculated by TB also seems to be less shifted in en- 
ergy with respect to the spin-down part as compared to 
the LAPW reference. This feature of DOS corresponds 
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FIG. 2. (Color online) Comparison of TB results assuming 
"atomic parameters" and second nearest neighbours (a,2nn) 
to LAPW using LSDA and GGA approximations for quanti- 
ties: (a) Spin magnetic moment on Co (thick lines) and Pt 
(thin lines) sites per formula unit with magnetization along 
the [001] axis (spin in the interstitial region of LAPW is 
not captured); (b) DOS per formula unit projected on spin- 
up/down states with magnetization along the [001] axis, the 
legend of plot (a) applies including the line style used for Pt 
and color coding; (c) MAE = E to t{M 110 ) - E to t{M 0Q1 ) per 
formula unit for a range of band fillings where the natural 
band filling is 19 valence electrons. 



to the observation of lower spin magnetic moment pre- 
dicted by TB both for Co and Pt at the Fermi level in 
Fig. [2]Ja) and motivates us to increase the Stoner param- 
eters both for Co and Pt atoms and add the third nearest 
neighbours to enhance the coupling between Co and Pt 
in the next section. 

Fig. [2^c) compares MAE calculated again using TB 
with "atomic parameters" and second nearest neigh- 
bours and LAPW with LSDA and GGA approxima- 
tions. Note that in the TB case MAE amounts to the 
difference of total energies for two magnetization direc- 
tions: MAE = Etot(Muo) — Etot(Mooi), whereas in case 
of LAPW we use the force theorem following the work 
of ShickPl (We use about 50000 and 200000 fc-points in 
LAPW and TB integration, respectively.) MAE is a more 
subtle quantity than the spin magnetic moment or DOS 
so we calculate it for a wide range of band filling (b.f.) us- 
ing the rigid band approximation and compare the trends 
of MAE(b.f.) rather than individual values. 

We observe a very good agreement of the LSDA and 
GGA data and a broad agreement of the TB and ab 
initio curves. We consider this a remarkable success of 
the TB model which employs parameters fitted for pure 
single-element metals. The amplitude of MAE oscilla- 
tions decreases towards the edges of the valence band so 
the relative error of the TB result at the Fermi energy 
becomes quite large. Both LAPW values MAElsda = 
1.26 meV/f.u. and MAE GGA = 0.85 meV/f.u. are in 
good agreement with the measured value MAE exp . = 
1.0 HieV/f.u.jS whereas our "atomic parameters" TB 
prediction based on straight forward transfer of Shi and 
Papaconstantopoulos' parameters is an order of magni- 
tude lower: MAE TB = 0.13 meV/f.u. 



V. PARAMETER REFINEMENT AND MAE 

The good agreement of key features in the spin mag- 
netic moment, DOS, and MAE calculated by LAPW and 
TB motivates us to improve our TB model. The short- 
comings of the initial simulations described in the previ- 
ous section provide useful guidance how to proceed. As 
concluded in the discussion of Fig. [2^b) , we add the third 
nearest neighbours to enhance the coupling between Co 
and Pt sites and we increase the on-site Stoner parame- 
ters both for Co and Pt atoms to match the larger spin 
polarization obtained by LAPW. 

The results are summarised in Fig. [3] which presents 
the same quantities as Fig. [2] computed in three dif- 
ferent ways: TB with "atomic parameters" and third 
nearest neighbours (a,3nn); TB with "atomic param- 
eters" enhanced by I Co ,d = 1.08Ig,* o ^f c = 0.134 Ry, 
Ipt.d = 0.15/co,d, and third nearest neighbours (a,I,3nn); 
LAPW with LSDA (same as in Fig. § GGA is not in- 
cluded to maintain legibility of the plots). 

Fig. [3ta) can be contrasted with its counterpart, 
Fig. [2ja)T The spin magnetic moment assuming the third 
nearest neighbours is again underestimated by TB us- 
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FIG. 3. (Color online) Comparison of TB results assuming 
third nearest neighbours and "atomic parameters" (a,3nn) 
or "atomic parameters" with enhanced Stoner parameter 
(a,I,3nn) to LAPW using LSDA approximation for quanti- 
ties: (a) Spin magnetic moment on Co (thick lines) and Pt 
(thin lines) sites with magnetization along the [001] axis; 
(b) DOS projected on spin-up/down states with magnetiza- 
tion along the [001] axis, the legend of plot (a) applies includ- 
ing the line style used for Pt and color coding; (c) MAE = 
Etot (Mno) — Etot (Moot ) for a range of band fillings where the 
natural band filling is 19 valence electrons. 



ing the default "atomic parameters" . However, we man- 
aged to tune the Stoner parameters Ipt,d and Ico.d to 
achieve quantitative agreement of the spin magnetic mo- 
ment on Pt with the LSDA values throughout the whole 
valence band. The spin magnetic moment on Co cannot 
be brought to a full quantitative agreement with LSDA 
by tuning only parameters Ipt t d and Ico.d an d remains 
slightly lower than the LSDA or GGA reference at the 
Fermi energy. We address this deficit in the final refine- 
ment of our TB model. 

DOS in Fig. [3^b) shows significant improvement over 
Fig. [2jb) so the summation to third nearest neighbours 
seems to be more suitable for modelling the coupling be- 
tween Co and Pt. Moreover, the enhancement of the 
Stoner parameter increased the mutual shift of spin- 
up/down DOS to show closer match with LSDA, as ex- 
pected. 

In Fig. 3^c) we focus on band fillings closer to the nat- 
ural band filling (b.f. = 19 valence electrons). Including 
the third nearest neighbours does not cause any signifi- 
cant change in the overall MAE dependence on the band 
filling. Slightly more pronounced change corresponds to 
increasing the Stoner parameters but the main features 
of MAE are still in agreement with LSDA and GGA data. 
Note that the MAE at the natural band filling becomes 
negative due to the enhanced exchange interaction. 

We can conclude that our first TB parameter refine- 
ment attempt presented in Fig. [3] leads to a better agree- 
ment with LAPW in a broad range of valence band ener- 
gies, however, increasing Ipt,d and Ico,d and adding third 
nearest neighbours does not reproduce accurately the net 
spin magnetic moment predicted by LAPW at the Fermi 
energy and the MAE at the natural band filling. 



As we mentioned in Sees. Ill and IV the unoccupied 
p-states were added to the Harrison TB model by Shi 
and Papaconstantopoulos only to produce more realistic 
warping of the d-band in single-element crystals close to 
Fermi energy. Therefore, the position of the unoccupied 
p-states in the CoPt compound is the next natural sub- 
ject to scrutiny. 

We note that Pt offers more room for variation of the p- 
state on-site energy as its s-state is much lower in energy 
than the s-state of Co: {£pt, P -£pt,d) / {£pt.s-£pt,d) = 2.2 
whereas {£co, P -£co,d)/{£co,s-£co,d) = 1-02. Therefore, 
we can bring £pt, P closer to £pt, s (and to ept.d further 
below) without changing the order of the Pt on-site en- 
ergies. Such shift should increase the hybridisation of the 
p-states on Pt with the exchange- split d-states increasing 
the deficient net spin magnetic moment. 

Fig. [4] shows an overview of TB and LAPW results 
analogous to Fig. [3] The new TB data are calculated 
using "atomic parameters" with enhanced Stoner pa- 
rameter, third nearest neighbours, and two examples of 
the shifted Pt on-site energies: {ept, P — £pt,d)/{£pt,s — 
£pt,d) = 1-61 labelled as (a,I,pl,3nn) and (£pt, P — 
£pt,d)/(£pt,s-£pt,d) = 1.26 labelled as (a,I,p2,3nn). The 
LSDA data are the same as in Figs. [2] and [3] 

As expected, the greater proximity of p and cZ-states 
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E-E p [eV] 




Band filling - natural b.f. [els./f.u.] 



FIG. 4. (Color online) Comparison of TB results assuming 
third nearest neighbours, "atomic parameters" with enhanced 
Stoner parameter, and on-site energies £pt, P — £pt,d = 0.69 Ry 
(a,I,pl,3nn) or £pt, p — ept,d = 0.54 Ry (a,I,p2,3nn) to LSDA 
approximation for quantities: (a) Spin magnetic moment on 
Co (thick lines) and Pt (thin lines) sites with magnetization 
along the [001] axis; (b) DOS projected on spin-up/down 
states with magnetization along the [001] axis, the legend 
of plot (a) applies including the line style used for Pt and 
color coding; (c) MAE = E t ot(M lw ) - E to t{M oi) for a range 
of band fillings where the natural band filling is 19 valence 
electrons. 



increases the spin polarization of Pt deeper in the va- 
lence band, however, the effects compensate at the Fermi 
energy giving excellent agreement with LSDA data as 
shown in Fig.[4|a). At the same time, the spin magnetic 
moment of Co at Fermi energy scales with the shift of 
on-site energy £pt, P due to hybridisation with Pt and fi- 
nally reaches the LAPW values when the shift mentioned 
above is in the range £pt, p — £pt,d ~ 0.54 — 0.69 Ry. 

The agreement of DOS calculated by LAPW and TB 
with enhanced Stoner parameters shown in Fig. |3jb) is 
satisfactory. We include Fig. Qb) to demonstrate that 
shifting £pt, P causes only minor deviation from DOS ob- 
tained by LSDA at the very bottom of the valence band, 
whereas it further improves the agreement closer to the 
Fermi energy. 

Fig. |2|c) shows the main result of our work. The gen- 
eral trend of the MAE dependence on the band filling is 
very robust against small variations of the TB input pa- 
rameters and is in broad agreement with MAE calculated 
by LAPW. On the other hand, the value of MAE for the 
natural band filling turns out to be very sensitive to the 
input parameters. Remarkably, the MAE for the natural 
band filling is in quantitative agreement with the LSDA, 
GGA, and experimental values MAE s» 1 meV/f.u. when 
the on-site energy of the p-states on Pt is in the range de- 
termined by comparing DOS and spin magnetic moment 
to LAPW: e Pt , p - e Pt4 « 0.54 - 0.69 Ry. 

We have also explored the sensitivity of the above 
quantities to the variation of the Slater-Koster param- 
eters. However, replacing the geometric average by an 
arithmetic average to obtain V a g' pt or increasing the rel- 
ative magnitude of the Slater-Koster parameters to en- 
hance the Co-Pt hybridisation does not change the spin 
magnetic moment or DOS in a transparent way that 
would improve our physical understanding of the elec- 
tronic structure or the overall agreement with the LAPW 
results. 



VI. SUMMARY 

We have carried out systematic modeling of electronic 
structure and relativistic magnetic characteristics of bulk 
CoPt LIq structure using TB and ah initio methods. An 
emprirical multiorbital TB model following the Harrison 
approach with a parametrisation devised by Shi and Pa- 
paconstantopoulos was applied. We extended the model 
by adding an atomic spin-orbit coupling term to the on- 
site Hamiltonian blocks in order to account for the mag- 
netocrystalline anisotropy. We have focused on the MAE 
as a function of the band filling so that we could com- 
pare general trends of TB and LAPW rather than singe 
values for the natural band filling which are available in 
literature. 

We started our calculations by checking the validity 
of the model assuming the "atomic parameters" (param- 
eters optimized for single-element fee crystals with geo- 
metric averaging of Slater-Koster parameters between Co 
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and Pt) by comparing the TB band structure, spin mag- 
netic moment, DOS, and MAE to corresponding LAPW 
results obtained using the Wien2K program package. 
The broad agreement observed throughout the valence 
band with slightly deficient TB net spin magnetic mo- 
ment and smaller mutual shift of the projected DOS 
stimulated further development of the parametrisation 
focusing on the Stoner parameters and on-site energies 
of the virtual p-states. 

We continued by adding the interaction to third near- 
est neighbours to enhance the hybridisation between 
magnetic and spin-orbit coupled sites and varied three 
TB parameters (Ipt,d> Ico,d> and £pt,p) which extend the 
original Harrison parametrisation and their values were 
likely to require corrections after the transfer to an or- 
dered binary alloy. We compared the new set of results 
to ab initio predictions again. In other works, this pro- 
cess is typically replaced by simultaneous fitting of all 
TB parameters to APW band structures and total ener- 
gies, however, in our work we seek better physical insight 
into the spin-orbit coupling phenomena and transferabil- 
ity between different structures rather than precise agree- 
ment of TB and ab initio results. 

We found that the net spin magnetic moment increases 
with small enhancement of the Stoner parameters and 
with shifting of the Pt p-state on-site energy towards the 
rf-states as expected. Remarkably, the MAE obtained 
for a narrow range of these parameters, where the net 
spin magnetic moment and DOS reached the best agree- 
ment with the ab initio predictions, is in quantitative 
agreement with ab initio results. Such success motivates 
future investigations of the transferability of the model 
in other compounds or multilayers. At the same time our 
TB model is well suited for incorporation of the equilib- 
rium Green's function framework to calculate relativis- 
tic magnetotransport phenomena in structures contain- 
ing the magnetic compounds or multilayers. 
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Appendix A: Spin-Orbit Hamiltonian 

The spin-orbit Hamiltonian term has to be written in 
the basis of cubic harmonics for which the Slater-Koster 
parameters are derived in literature! 21 * 23 ! The cubic har- 
monics can be written in terms of the angular momentum 
cigenstates for the s, p, and d orbitals: 
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P.r 




Pv 


- — (Y- 1 


+ 1? 


Pz 


= *?, 




d X y 






dyz 




+ 1? 


d>xz 


= 7i0r' 


-Yi 


d>x 2 — y 2 


= 


+ Y 2 2 


d 3z 2_ r 2 


- Y° 





(Al) 

where (n\l,m) = Yj m (0, L z \l,m) — m\l,m), and 
L 2 \l,m) — 1(1 + l)\l,m}. In atomic units the angular 
and spin moments are dimensionless and the coefficient 
A aQ in Eq. ^ is measured in Rydbergs. With cubic har- 
monics ordered as in Eq. (Al) we obtain the following 



Hamiltonian contribution for the p orbitals: 
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and for the d orbitals: 
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(A3) 



where the first and second diagonal blocks correspond to 
spin up and down eigenstates: S z \l Sl m s ) ~ ±l/2|^ s , m s ), 
respectively. The total Hso matrix is added to the on- 
site Hamiltonian terms, whereas the inter-site matrix el- 
ements remain unchanged. The size of X aa depends on 
the type of atom. We neglect Xco,p and Xpt, p as they are 
at least an order of magnitude smaller than Xco,d and 
Xpt,d- 

Appendix B: Spin-Orbit Hamiltonian 

As mentioned in Sec. [TT] the hopping integrals 
Ebp,aa{pn) can be written in terms of Slater-Koster^ 



parameters VV™ (neglecting the site indices) for the s 
and p orbitals: 

E — V 

Es,x IVspai 

Ex,x — I Vppa ~t~ (1 ^ )^pp7ri 
Ex.y iTYlVppa linVpp^^ 

E x , z = InVppv - InVppn, (Bl) 
and for the s, p, and d orbitals: 



E s .xy = VSlmVsda, 

E StX 2_ y 2 = ~y^ 2 ~ m 2 )V sd a, 
E s , 3z 2_ r 2 = [n 2 ~^(l 2 + m 2 )]V sd<J , 

E x , X y = V3l 2 mV pdt7 + m(l - 2l 2 )V pd7T , 
E XiVZ = VsimnVpdo- - 2lmnV pdlTl 
E x>zx = V3l 2 nV pdrT + n(l - 2l 2 )V pd „, 

Vs 

E XlX 2- v 2 = -t^IQ 2 ~ m 2 )Vp da + l(l - I 2 + m 2 )V pd7r , 
Ey tX 2_ y 2 = —m(l 2 - m 2 )V pda - - m(l + I 2 - m 2 )V pd7r , 
E z , x 2_ y 2 = ^~n(l 2 - m 2 )Vpda - n(l 2 - m 2 )V pd7: , 
E x ,3z2~ r 2 =l[n 2 - ^{l 2 + m 2 )}V p da - V3ln 2 V pd7T , 
E y , 3z 2^ r 2 = m[n 2 - ~(J 2 + m 2 )]V pda - V3mn 2 V pd ^, 

E z ,z z 2- r 2 = n[n 2 - * (Z 2 + m 2 )]V pda + V3n(l 2 + m 2 )V pd7r , 

Exy.xy = 3l 2 m 2 V dda + {I 2 + rn 2 - 4l 2 m 2 )V dd7I + (n 2 + l 2 m 2 )V ddS , 
E xy ,yz = 3lm 2 nV ddl7 + ln(l - Am 2 )V ddn + ln(m 2 - l)V dd s, 
E xy ,zx = il 2 mnV ddcJ + mn(l - 4l 2 )V ddn + mn(l 2 - l)V ddS , 
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E 



E 2 2 — 

^yz.x^ —y 



E 2 2 

^zx.x^ —y 



Ex 2 — y 2 ,x 2 — y 2 



3 1 

>xy,x 2 -y 2 = ^m{l 2 - m 2 )V dd(J + 2/m(m 2 - l 2 )V dd7T + -lm(l 2 - m 2 )V dd&1 

3 1 

-mn{l 2 - m 2 )V dda - mn[l + 2(l 2 - m 2 )]V dd7[ + mn[l + -(I 2 - m 2 



2 V " •■"')WddS, 

^nl(l 2 - m 2 )V dda + nl[l - 2{l 2 - m 2 )]V dd7r - nl[\ ~\{l 2 ~ m 2 )]V dd5 , 



E X y,3z 2 -r 2 = VO 
Ey Z ,3z 2 -r 2 = VO 
E Z x,3z 2 -r 2 = Vo 



lm{n 2 ~\{l 2 + m 2 ))V dda - 2lmn 2 V dd7I + -Zm(l + n 2 )^ 
mn(n 2 - ^(Z 2 + m 2 ))^^ + rnn(7 2 + to 2 - ?i 2 ) VddTr - ^mn(l 2 + m 2 )V dd s 



Ex 2 —y 2 ,3z 2 — r 2 



ln (n 2 - \ ( l2 + m 2 ))V dda + ln(l 2 + m 2 - n 2 )V ddn - ^ln(l 2 + m 2 )V ddS 
= \{l 2 - m 2 ) 2 V dda + [I 2 +m 2 - (I 2 - m 2 ) 2 ]V dd „ + [n 2 + \{l 2 - m 2 ) 2 ]V ddS , 
V3 r ' 



E3Z 2 -r 2 ,3z 2 -r 2 = [n 2 - - (I 2 



^(l 2 - m 2 )[n 2 -\{l 2 + m 2 )\V dda + n 2 {m 2 - l 2 )V dd7T + + n 2 )(l 2 - m 2 )V t 
(I 2 + m 2 )] 2 V dda + 3n 2 (l 2 + m 2 )V ddv + \{l 2 + m 2 ) 2 V ddS , 



(B2) 



where I, m, and n are the directional cosines of the inter- denote the cubic harmonics written in terms of spherical 
site position vector p n and the labels s, x, y, z, xy, etc. harmonics in Eq. (All. 
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